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We study a O-tt dc superconducting quantum interference device (SQUID) with asymmetric in¬ 
ductances and critical currents of the two Josephson junctions (JJs). By considering such a dc 
SQUID as a black box with two terminals, we calculate its effective current-phase relation Isijp) 
and the Josephson energy where if is the Josephson phase across the terminals. We show that 

there is a domain of parameters where the black box has the properties of a JJ with degenerate 
ground state phases ip = ±ip. The ip domain is rather large, so one can easily construct a ip JJ 
experimentally. We derive the current phase relation and show that it can be tuned in situ by 
applying an external magnetic flux resulting in a continuous transition between the systems with 
static solutions ip = ±(p, ip = (po (ipo 7 ^ 0, tt) and even ip = poAnp. The dependence of po on applied 
magnetic flux is not An (one flux quantum) periodic. 

PACS numbers: 74.50.-|-r, 85.25.Cp 
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I. INTRODUCTION 

One of the most important characteristics of a Joseph¬ 
son junction (JJ) is its current-phase relation (CPR), i.e., 
the relation between the supercurrent Ig flowing through 
the junction and the Josephson phase (p across it. The 
CPR, aka the first Josephson relation, plays a fundamen¬ 
tal role and is responsible for almost all properties of the 
JJ. The CPR usually depends on the microscopic physics 
of the Josephson barrier. The CPR is usually a 27r peri¬ 
odic function of (p and, in the most simple case, is given 
by Is = Ic sin (p. A review of different types of CPRs can 
be found elsewhere^. 

The quantity directly related to the CPR is the Joseph¬ 
son energy profile U{(p). It is defined so that Is{(f>) = 
dU{(p)/d(p. For a sinusoidal CPR U{cp) = Ej{l — cost/)), 
where Ej = ^oIc/2Tr is the Josephson energy, and 
<i)o Ri 2.068 fWb is the magnetic flux quantum. U{(p) 
is also a 27r periodic function of (p. 

The most interesting JJs both in terms of applications 
and fundamental studies are those having a non-trivial 
Josephson energy profile U{(p). For example, tt JJs that 
received a lot of attention^^^^ have negative critical cur¬ 
rent, which results in the sign change of U{(p). As a 
result, the energy minimum (and the ground state) cor¬ 
responds to (/) = TT, where the usual JJ (called 0 JJ in 
this context) has a maximum. Further, there was a theo¬ 
retical predictions^^d4 obtain the so-called 

ipo JJ, i.e., the one having a single U{<p) minimum (with 
each 27r period) situated at (p = ipo 0, tt. More recently, 

JJs having a periodic double-well potential U{<p) and, 
therefore, a degenerate ground state (two U{(p) minima) 
with the phases (p = Eip were investigated^®^^®. Upon 
application of magnetic field one can also obtain a U{(p) 
without reflection symmetry. 

Theoretically, it was predicted^®’^° that a JJ made of 
d-wave superconductors in a specific range of parame¬ 
ters (orientation angle, temperature, etc.) can posses a 


ip ground state. A degenerate ground state was obtained 
from measurements of the CPR of d-wave based nano 
JJs^^. Later on, some other indications of a ip state, such 
as an anomalous temperature dependence of the critical 
current, were observed^^ also on nano JJs. The faceting 
along longer grain-boundary JJs based on d-wave super¬ 
conductors results in an effective (facet-averaged) ip JJ^^. 
In the latter case one expects non-quantized splintered 
vortices^"^’^®, which were observed^® using superconduct¬ 
ing quantum interference device (SQUID) microscopy. 
However, in all cases no state manipulation or readout 
were demonstrated, probably because of high damping 
and poor control over JJ properties. 

Recently, we have suggested^® and successfully 
demonstrated^® a ip based on conventional low-Tc su¬ 
perconductors with tailored ferromagnetic barrier^^’^®. 
This junction has a degenerate ground state phase Eip, 
i.e. its Josephson energy profile looks like a 27r-periodic 
double well potential (in the absence of bias current). 
The two ground states can be used to store information^^. 
The unusual physics of ip JJs was discussed in several 
works^®’^®. However, the ip JJs constructed so far are 
rather large. Making them smaller (shorter) requires very 
exact control on parameters such as critical current den¬ 
sities and the lengths of the 0 and tt regions. 

In this paper, we propose an effective ip Ii based on 
an asymmetric O-tt dc SQUID, i.e., a dc SQUID with one 
0 and one tt JJ, with finite inductance and asymmetric 
critical currents of the 0 and tt JJ. This not only has 
advantages over the previous proposals^®’^®’^® in terms 
of geometrical dimensions, margins, and the size of the ip 
domain in parameter space, but also shows other unique 
features. For example, it can be operated not only as a 
ip JJ, but also as a ipq JJ^® or as a combination of both, 
i.e., as a ip^Eip JJ. This goes much beyond the trivial use 
of a dc SQUID as the substitution for a single JJ with 
magnetic-field-tunable critical current. 

The paper is organized as follows. In sec. H we intro- 
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duce the considered model system and the equations de¬ 
scribing it. In the main sec. Ill we present our numerical 
results. Sec. IV summarizes our findings. In appendix A 
we consider the vicinity of the 0 •<->■ (p transition in param¬ 
eter space and derive many results analytically similar to 
the Ginzburg-Landau approach. 

II. MODEL 

The SQUID circuit is shown in the inset of Fig. I. The 
applied bias current / splits into two branches with induc¬ 
tances Li and L2 and Josephson junctions with critical 
currents Id and Ic2- 

For the sake of simplicity, we will derive everything in 
normalized units. The current will be normalized to the 
largest (by absolute value) critical current of the two JJs. 
Without loss of generality we assume that Id > 0 and 
\^c2\ ^ Id • Then we define (x — Ic 2 /Ic\ the normalized 
critical current of the second JJ, with |Qf| < 1. Thus, 
a > 0 corresponds to a 0-0 SQUID, while a < 0 to a O-tt 
SQUID. Although, our main focus is on a O-tt SQUID 
(a < 0) the results for a 0-0 SQUID will be included 
automatically as well. Further, we introduce normalized 
inductances 

2-KldLi „ 27r/ciL2 , . 

A = A = (1) 

Note that the definition of /32 uses Id , so that by chang¬ 
ing Pi and P 2 one can see the effect of inductances only, 
while changing a one can see the effect of critical cur¬ 
rent asymmetry only. These definitions are related to the 
conventional^ Pl = 2 /ciTe/^'o, as ttPl = Pi + P2 = Ps, 
where Ls = Li + L 2 . 

The total phase ip across the SQUID, see the inset of 
Fig. 1, can be expressed in two ways 

ip = (pi+ Pisincpi + ri(pe\ ( 2 a) 

Ip = (p2 + aP2Sm(p2 - r2(pe, ( 2 b) 

where ri -|-r 2 = I, and ri and r 2 are the ratios, in which 
the externally applied normalized flux/phase (pe = 27r/ = 
27r<I>e/‘ho (/ is the normalized flux, aka frustration) is di¬ 
vided between the two branches. If the external magnetic 
field is applied by a coil, ri^ 2 <Pe = 27rMi^2do/$o, where Iq 
is the current in the coil creating the magnetic field and 
Mi ^2 is the mutual inductance between this coil and the 
left and the right arms of the SQUID, respectively. For 
the sake of simplicity, we neglect the mutual inductance 
between Li and L2. 

The sum of the currents in both branches is given by 
7 = sin -I-a sin (/) 2 , (3) 

where "f = I/Id is the normalized bias current. 

Since we are interested in a geometrically small system 
which should not have too many internal states, we focus 
on the case of small, but finite, inductances, i.e., 0 < < 

1 and 0 < P 2 < I- 



FIG. 1. (Color online) The energy U{ip) of the system given 
by Eq. (5) for different values of the asymmetry parameter 
a and for Pi = P 2 = 0.7 (qc ~ —0.417). Inset shows the 
schematics of the circuit considered. 

III. RESULTS 

The easiest way to solve Eqs. (2) and (3) is to cal¬ 
culate 01 and 02 from Eqs. (2) for given ip. Note that 
for Pi,P 2 ,\a\ < 1 each of Eqs. (2) has a unique solu¬ 
tion 01 or 02 for given ip. It is convenient to define a 
universal function <p{v,p), which is a solution 0 of the 
equation v = <p + p sin( 0 ) for given argument v and pa¬ 
rameter \p\ < 1 . Then (piiip) = 0(0 — ripe, Pi) and 
02 ( 0 ) = 0(0 + '>’ 2 pe,OiP 2 )- The function p{v,p) has an 
obvious, but useful, property: 

0(u-b 7 rn,p) = 7 rn-I-0[n, (-I)”p], (4) 

where n is an integer. 

A. Current-phase relation and Josephson energy 

To determine the CPR 7 ( 0 ) we calculate 0i(0) and 
02(0) as mentioned above and then use Eq. (3) to obtain 
7 . We do not show 7 ( 0 ) plots here, but show Josephson 
energy U{ip) plots instead. 

The total energy of the system is given by 

U{ip) = Uj{ip)+ Ul{iP), ( 5 ) 

where 


Uj{ip) = [1 - cos01 (0)] -b a;[I - cos02(0)], (6) 
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FIG. 2. (Color online) The domain of the ip state (pink/gray) 
as a function of parameters a, j3i and ^ 2 - 

is the Josephson energy of both JJs, and 

UlW = Y sin^ sin^ (7) 

is the magnetic field energy stored in the inductors. By 
direct substitution, one can see that = 'y{' 4 ’) (here 

and below the prime denotes 8 /dip by default), like for 
any JJ. 

Consider the case of zero applied magnetic flux = 0. 
Several examples of U{ip) are presented in Fig. 1. At 
positive a (conventional 0-0 SQUID) the energy profile 
resembles the usual UconvW = 1 ~ cos(^) profile of a 
conventional single JJ with the ground state at "0 = 0. 
Deviations from 17conv(V') due to finite inductances 
and make U{ip) sharper than C/conv('*/') near the maxima 
and more shallow than C/conv ('*/') near the minima. As 
a decreases down to 0, the height of U{ip) decreases as 
2(1-I-a) while the shape of U{ip) almost does not change. 
As a becomes negative (O-tt SQUID) the U{ip) becomes 
more flat near ip = 0. For a below some critical value 
ttc the energy profile develops two minima, as can be 
seen in Fig. 1. The system has a degenerate ground state 
Ip = ±ip, j.e., forms a (p jjis,16,18-20,23,25,31^ 

The transition to the degenerate ground state takes 
place at the value of a = ac, for which U"(Q) = 7^(0) = 0. 
This allows to calculate 


I + 81 + 132 1 -I- /3e 1 + TT/Si 

This is one of the central results of the paper. 

The domain of the p state is shown in Fig. 2. From the 
sides it is limited by our choice of parameters 0 < /3i,2 < 1 
and from the bottom by a > —1. From the top it is 
limited by the surface ac{l3i, 82 ) given by Eq. (8). The p 
domain has the maximum size (height) along the a axis 
equal to —1 < a < —1/3 for /3i_2 = 1- For 81^2 —>■ 0 the 


height of the p domain along the a axis vanishes linearly 
oc /3s- We also note that in the case of /3 e = 0, used 
by many authors because it is solvable analytically, one 
never obtains a 7; JJ. 

An important practical difference between the 0- 
TT dc SQUID considered here and a O-tt JJ of finite 
length^®’^®’^^’^^ is that our system has a rather large p do¬ 
main where the theory presented here works. Even in the 
case of /3i; —>■ 0 the range of a = — 1.. .ac, correspond¬ 
ing to the p domain, shrinks as Aa = ac — (—1) = /3e, 
le., linearly. Instead, for O-tt JJs the p domain shrinks 
asi5.3i g ~ ^ where 5 = (Lq — L,r)/2 is the devia¬ 

tion of the 0 and tt facet length from the average length 
L = [Lq -f Lt^)/2. This imposes additional requirements 
on fabrication accuracy or one should move to a region of 
the phase diagram where the theory works only qualita¬ 
tively, as in the first experiments^®’^^. In other works®^’®® 
some alternative techniques to enlarge p domains were 
suggested. However, the relative volume of the p do¬ 
mains in parameter space is still much smaller than in 
the present work. 

Let us now consider U{ip) at finite applied magnetic 
flux (pc ^ 0. First, we point out that Eqs. (2) and (3) 
have several important symmetry properties, when the 
applied flux changes by a half-integer number of flux- 
quanta $0 

(p1^^ = (pe + mr; (9a) 

(9b) 

= </>i; (9c) 

(p 2 ™ = (p 2 + mr- (9d) 

^new ^ (9e) 

(9f) 

Upon such transformation the current does not change, 

ie., 

7(^, a) (10) 

while the Josephson energy, which is still 2tt periodic in 
Ip, can be expressed in terms of new variables as 

U{iP, a) = U{iP’^^'^, a“"') -f «“"'[(-!)’" - 1]. (II) 

Thus, by applying a half-integer / we can turn a 0-0 
SQUID with a> \ac\ (effective 0 JJ) into a O-tt SQUID 
(effective p JJ) with = —a < ac and vice versa. 
An even more interesting point is that ip shifts by an 
amount rimr, which is, in general, not a multiple of 27r. 
Since after transformation (9) the value of 7 does not 
change (it depends only on {(pi,(p 2 ) mod 27r), the CPR 
"/(ip) shifts along the ip axis by the amount rim:. This 
is another central result of the paper. 

Examples of U{ip) calculated for different values of ap¬ 
plied magnetic flux are shown in Fig. 3. For a > \ac\, 
see Fig. 3(a), U{ip) at / = 0 has a single minimum at 
Ip = 0. At / = 0.5 it transforms into a double well po¬ 
tential. It is the same potential as one gets for a = —0.5 
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FIG. 3. (Color online) Josephson energy !/('(/)) for — 0.4, 
/?2 = 0.6, ri = 0.4 and different / = 0... 1 specified next to 
each curve, (a) a = 0.7 > \oic\ = 0.5; (b) a = 0.4 < |ac| = 

0.5. 


and / = 0, but shifted by (centered at) ip = riw = 0.47r 
and lifted by 2a = —2a"®"' = 1.4, see Eqs. (9) and (11). 
In essence, at this / one obtains a (po ^ ‘P JJ- For / = 1 
the energy U{ip) transforms again to the same profile as 
at / = 0, but shifted by (centered at) ip = ri27r = O.Stt, 
see Eqs. (9) and (11). For 0 < a < |ac|, see Fig. 3(b), 
we again start at / = 0 from having a single minimum 
of U{ip) ai Ip = 0. At / = 0.5 we again obtain the U{ip) 
profile as the O-tt SQUID with a"®"' = —a, but at / = 0, 
shifted hy Ip = tity. However, in this case it is not a 
double-well as —a ^ ac- At / = 1 we obtain the original 
profile (as for / = 0), but shifted by ri27r, see Eqs. (9) 
and (11). 



FIG. 4. The ground state phase ±(p(Q:). The horizontal 
dashed line shows the unstable static solution ip = 0 in the 
region a < a^. Thin line shows the approximation given by 
Eq. (All). 


rather simple transcendental equation, 

sin((^i) -|-asin[(/)i -|-sin((^i)] = 0, 

which we solve to find pi. Then the value of p 2 , if neces¬ 
sary, can be calculated from Eq. (12). The ground state 
phase tp is obtained using one of the Eqs. (2), with /3i 
and P 2 corresponding to the real circuit. 

An example of p{a) at fixed ( 3 i ^2 is shown in Fig. 4. 
In essence this is p along a vertical line crossing a p 
domain in Fig. 2. The ground state phase is zero as a 
decreases from 1 down to the bifurcation point ac- After 
the bifurcation point (a < ac) the zero solution becomes 
unstable, as indicated by the dashed line. Instead two 
degenerate stable solutions appear. 

The phase p corresponding to the degenerate state de¬ 
parts from zero as described by Eq. (All). At a —>■ —1 
the ground state phase tends to its maximum value 
which is given by 


_ TT ^ Pi- P2 

< 7 >max- 2 +^2’ 

where y is a solution of the equation 


(13) 


B. Ground state phase 

The ground state (/ = 0, 7 = 0) phase p in the 
p domain a < ac can be calculated numerically from 

^{p) = 0 . 

To simplify and accelerate this computation we use the 
following procedure. Since for 7 = 0 the phases pi and 
p 2 (but not p) depend only^^ on /3s, for calculation of 
pi and p 2 , without loosing generality, we assume that 
Pi = /3e and P 2 = 0. Then Eqs. (2) collapse to 

02 = <(>1+/3Esin((^i). (12) 

Substituting this into Eq. (3) with 7 = 0 we arrive at a 


2y = PsCOs{y). (14) 

The possible range of y is from j/ = 0 for /3i = (32 = 0 
to y = y*, where y* « 0.739 is a solution of the equation 
y* = cos(y*). The corresponding range of '0max is then 
from tt/2 — y* (reached for Pi —>• 0, (32 —>■ 1) to ty/2 + y* 
(reached for /3i —>• 1 and P 2 —>• 0). Note that since y 
is a function of /3 e only, it follows from (13) that the 
phase i/?max is antisymmetric with respect to the diagonal 
direction Pi = P 2 - In particular, pmax = § at /3i = 02 
(symmetric system). It is interesting that even for a tt JJ, 
which is weaker than the 0 JJ, one can obtain a ground 
state phase \p\ > tt — a situation, which is not possible 
in a continuous O-tt JJ studied earlier^^J®. 
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FIG. 5. The phase diagram of O-tt dc SQUID (ground state phase (/p(/Ii,/ 32) for (a) a = —0.6 and (b) a = —0.8. The dashed 
line shows the boundary between trivial '0 = 0 and tp = ground states given by expression (15). Continuous lines are the 
lines of the constant ground state phase tp. Its value is given next to each line. 


It is also interesting to plot the ground state phase 
on the (/3 i,/ 32) plane for fixed a, i.e., in essence, in the 
horizontal plane crossing the ip domain in Fig. 2 at fixed 
a. In such a plane, with the help of Eq. (8), the ip domain 
is given by 

02 >-1 — 01 for a < 0. (15) 

a 

Thus, the boundary is just a straight line, see Fig. 5. 
Below this line the ground state phase is 0 = 0, while 
in the filled area above this line the ground state phase 
is 0 = ±ip. Note, that the boundary given by Eq. (15) 
shifts towards the origin as a decreases. At a = — 1 the 
domain of a 0 state vanishes completely. At a < — 1 
(a case not considered here) a small ip = -k domain will 
appear close to the origin and will grow as a decreases 
further. 


C. Persistent current 

If the ground state of the system (/ = 7 = 0) is the ±ip 
state, one has a persistent current circulating clockwise 
or counterclockwise around the SQUID. From Eq. (3), its 
value is given by 

/circ = sin[0i(±(/j)] = asin[02(±</5)]- (16) 

Since in the ground state the phases 0i and 02 depend 
only^"^ on 0s, so does Icirc- This means that the value ip 
of the ground state phase, which depends on 0i — 02 and 
0 E, can be chosen independently from the value of the 
persistent current /circ, which depends only on 0s. 



FIG. 6. Girculating current as a function of a. Tilted dashed 
line shows the value of the critical current of the a-junction. 
The horizontal dashed-dotted line shows the maximum value 
of persistent current reached at a = — 1 and calculated using 
Eq. (17). Vertical dashed-dotted line shows the value 
calculated using Eq. (18). 

For fixed 0s the value of persistent current grows oc 
\Jac — OL near cuc and reaches the maximum value at a = 
— 1, see Fig. 6. This maximum value is equal to sin 02, 
where 02 is a solution of the following transcendental 
equation 

202 — 0s sin02 = TT. (17) 

In the limit 0s —>■ 0, from Eq. (17) 02 —>■ 7r/2 and circu¬ 
lating current —^ 1. 

Another interesting question is: does the value of circu¬ 
lating current ever equal to the maximum possible value 
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|a| — the critical current of the weaker JJ? It turns out 
that for fixed /3 e this happens for a = which is a 

solution of the following transcendental equation 

o-Pt. + arcsin(a) + ^ = 0, (18) 

see also Fig. 6. To proof this we take the state with 
02 = 7’'/2, then sin 02 = 1- From Eq. (3) sin0i = —a > 
0, and therefore 0i = arcsin(—a) (exactly this root!). 
By substituting this into Eqs. (2) and, without loosing 
generality assuming /3i = 0, 02 = we arrive at (18). 


D. Self-generated flux 

The self-generated flux in the loop (not including ex¬ 
ternal flux 0 e) is given by 

$ 

27r— =01 sin0i('0)-a02sin02('0) = 02(V^)-0i(0)-0e- 

(19) 

Similar to the circulating current, the value of sponta¬ 
neous flux in the absence of the bias (7 = 0 ) depends 
only^"^ on 0 s, which is not obvious at all from the first 
part of Eq. (19), but apparent from its last part. This 
allows to write 

$ 

27r— = /3s sin 01 = -a/3s sin 02 = 0s-fcirc- (20) 
^0 

Therefore, spontaneous flux behaves similar to the circu¬ 
lating current. The maximum value of the flux is also 
reached at a —>• — 1 . 
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FIG. 7. The amplitudes of the spontaneous circulating cur¬ 
rent /circ( 0 E) and spontaneous flux $( 0 e)/$o for fixed a = 
-0.7. 


a maximum, where Ichc = |ct| — the maximum possible 
value in our SQUID. It happens at 0 e = (2arcsin(a) -|- 
7 r)/(—2a), which was obtained by inverting Eq. (18). 


E. Critical currents 


In general, our system has four critical currents in the 
ip domain and two critical currents outside of it. These 
critical currents correspond to the escape of the phase 
from the left {—p) or the right (-ftp) wells of the double¬ 
well potential, see Fig. 1, for different directions of the 
bias current 7 . The critical current corresponds to the 
maximum of the CPR 7 ( 0 ). 




FIG. 8. Numerically calculated dependences 7 c±(Q!) for (a) 
01 = 02 = 0.7 and (b) 0 i = 0.1 and 02 = 0.9. 


It is instructive to replot Idrc and d> as a function of 
0E at fixed a, see Fig. 7. One can see that at small 
inductance 0 e < (1 -I- a)/{—a) [inverted Eq. (8)] the 
system is in the zero state. At larger inductances the 
system enters into ±p state and the spontaneous flux and 
current increase. However, $(0x;) grows monotonously, 
reaching ^ $ 0/5 at 0i; = 2, while the /circ(0E) exhibit 


The numerically calculated dependence of 7 c (a) for 
given 01 and 02 is shown in Fig. 8 . One can see several 
key points on the 7 c±(a) dependence. First, 7 c-i-( 0 ) = 1. 
Second, for 0i = 02 , 7 c-(—1) = 7 c+(—1), ie., the two 
dependences converge at a = —1, see Fig. 8 (a). In the 
case 01 < 02 , see Fig. 8 (b), the dependences cross at 
some — I < Ox < ttc- 
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F. Critical current as a function of applied flux 


The presence of the degenerate ground state and two 
critical currents also manifests itself in the dependence of 
the maximum supercurrent on magnetic field 7c(/)- In 
the general case this dependence can be calculated only 
numerically. Several examples of 7c(/) dependences are 
presented in Fig. 9. For the parameters in the (p domain 
(a < Qfc) one sees the domains corresponding to the dif¬ 
ferent flux states overlapping in the vicinity of integer 
/. In this overlapping region one observes in total four 
critical currents corresponding to the escape of the phase 
from different energy minima (wells) of the Josephson 
potential in two different directions. Note the striking 
similarity of these curves with those for JJs based on 
continuous O-tt JJs^®. 

As a increases and approaches ac, see Fig. 9, the 
bistability region is shrinking. Exactly at a = ac the 
overlap near integer / disappears, which corresponds to 
the disappearance of two distinct ground states ^|J = 
±(/ 5 . However, for a small range of a > ac, one ob¬ 
serves a small triangular-like bistability region where the 
branches meet, see Fig. 9(d). Further investigation of this 
small domain is outside the scope of this paper. We note 
that the intersecting domains, similar to those shown in 
Fig. 9 were calculated long ago^^’^®. At that time, how¬ 
ever, TT junctions were unknown, so only a 0-0 SQUID, 
where the domains intersect near half-integer /, was con¬ 
sidered. However some key results can be adopted to our 
case easily. 

In particular, we can find the positions of key points 
±A and zLX where d^c{f)/df = 0. We note^®’^® that 
the phases (j)i = ±7r/2 and 412 = ± 7 r /2 (in any combi¬ 
nation) satisfy the equation 7 ' = 0 for any asymmetry 
and applied flux. It turns out^® that each of these four 
combinations corresponds to an extremum of the 7 (/) de¬ 
pendence, ie., to points ±A and ±X in Fig. 9. To find 
the values of the external flux 27r/ = (j)e, correspond¬ 
ing to these points, we substitute the above phases into 
Eq. (2) and obtain 


= ± 


TT 

12 


Pi 


T 


TT 

12 


a/32 


( 21 ) 


The corresponding values of the critical current are ob¬ 
tained from Eq. (3) 


yextr = ±1 ± (22) 

In Eqs. (21) and (22) both ± signs are independent, 
providing four combinations in total. The trivial con¬ 
sequence of Eq. (22) is that for a SQUID with symmetric 
critical currents (|a| = I) the points ±X are situated at 
the horizontal axis ( 7 c = 0 ). 

The bistability region around integer / can be used to 
store one bit of information in the iip states as demon¬ 
strated recently^^. In some sense it is similar to the ear¬ 
lier proposals^^’^® to use n = 0 and n = 1 states of the 
SQUID biased to / ~ 1/2. However in our case, the flux 
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FIG. 9. The dependence of the critical cnrrent jc of the device 
on the normalized applied magnetic field (frustration) / for 
Pi = P 2 = 0.7 {nPh ~ 1.4, oic ~ —0.417) and different a. 


bias is not needed. Asymmetry also provides different 
critical currents at / = 0 , which simplifies readout. 

Finally, the practical question is: can one detect ±^ 3 , 
ifQ and ipa ± ip states in experiment by measuring the 
7c(/) dependence? For the state the answer is given 
in Fig. 9. One should observe intersection of branches 
and four critical currents. For the ± p state the sit¬ 
uation is similar. However, since this state appears at 
half-integer / the whole 7c(/) curve is shifted so that 
the bistability regions are situated around half-integer /. 
Finally, in the po state, which in our system appears at 
finite field only, one has only two critical currents and 
the junction looks just like a conventional one, although 
the 7c(/) dependence is unusual (periodic), but never 
multi-state. To prove the pq state one has to do a phase 















sensitive experiment, e.g., putting our black box in a su¬ 
perconducting loop. 


IV. SUMMARY 

We have shown that an asymmetric O-tt SQUID can 
be used as an effective </? JJ with magnetic field tunable 
current-phase relation 7('0) and, accordingly, a Joseph- 
son energy U{'iji). The critical value Oc of the critical 
current asymmetry parameter a required to obtain the 
degenerate ground state depends on the sum Li +L 2 
of inductances in two branches of the SQUID. Upon ap¬ 
plying an integer number of flux quanta f = n, the phase 
'i/' across the structure advances by an amount ri ■ 27rn, 
where r\ is the fraction of external flux induced in the left 
branch of the SQUID. Since, in general, ri is an arbitrary 
number depending on design, the phase shift ri • 27rn is 
not a multiple of 27r. By applying a half-integer number 
of flux quanta / = n-|-l/2toa 0-0 SQUID with a < |ac|, 
one can turn the effective 0-JJ into & ipa = riirn JJ. If 
a > \oic\, then one can turn the effective 0-JJ into a JJ 
with ground states The dependence of the critical 

current on magnetic flux clearly shows bistability regions 
typical for ip JJ^^J®. 

In terms of designing a practical device (bistable JJ) 
the target parameters can be, e.g., /3i = ,02 = 0.4... 0.7 
and a ~ —0.7... — 0.8 to stay well inside the p domain 
in Fig. 2. This will provide very large operation margins. 
Note that a finite inductance is essential to obtain a p 
domain of finite size. In the limit /3 e 0 the i^-domain 

shrinks to a point. 
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Appendix A: Solutions close to the bifurcation point 

a « Oc 

The JJ proposed here can also be used as a qubit 
at / = 0 , when the barrier separating the —p and +p 
states is small. This is the case when a is only slightly 
smaller than ac- In this limit the expressions for many 
important quantities can be obtained analytically. We 
are especially interested in the situation near the bottom 
of the energy profile U{gp), he., for small values of ip since 
the important physics (formation of ±(/3 state, escape, 
macroscopic quantum tunneling) takes place there. 

For a = Oc — e (0 < e <C 1) the Josephson energy of 
the system can be expanded like in the Ginzburg-Landau 


theory as 

Ugi,{iP) K. aip"^ ^-hip"^, |^/>| < 1. (Al) 

where the coefficients a{a) and h{a) have to be deter¬ 
mined from our model given by Eqs. (2) and (3), namely, 
from 

[/"(O) = 7'(0) = 2a; (A2) 

U""{Q) = 7'"(0) = 246, (A3) 

From Eq. (3) 

^'{ip) = cos{pi)(p[ + acos{p 2 )(p' 2 - (A4) 


The derivatives p[ 2 can be calculated by differentiating 
Eqs. (2) with respect to ip: 


1 

1 + ^1 cos((/)i) ’ 


1 

1 -I- j32acos{p2)' 


(A5) 


By substituting <p'i 2 from Eq. (A5) into Eq. (A4) we ob¬ 
tain 


7'(V’) 


cos( 0 i) 
l-b/3i cos((/)i) 


acos{p 2 ) 

1 -I- j32acos{p2)' 


(A 6 ) 


According to Eqs. {2), ip = 0 corresponds to (pi = p 2 = 
0, so from Eqs. (A2) and (A 6 ) we obtain the explicit 
expression for a(pt): 


a{a) 


1 + Pi 


a 

1 -I- 0/32 


(A7) 


which is negative for a < ac- Near ac the leading term 
is 


a(e) 


1 1 

"2 02(1+/3i)2""“'"- 


(AS) 


Similarly, differentiating j'{ip) in Eq. (A 6 ) two addi¬ 
tional times and using Eq. (A5) after each differentiation, 
we obtain a.t ip = (pi = p 2 = 0 from Eq. (A3) 


Ha) = id 


1 


a 


24 L(l + ^i)^ (l+o^ 2 )^J 

Near Oc the main term of 6 ( 0 ) is a constant 


(A9) 


./ x _ 1 fe(/3| + 3/3s -I- 3) ^ 1_ Q!^ -I-1 ^ 

™“24 (l + /3i)4 - 24 (-Oc)3(1 + /3i)4 

(AlO) 

Now we can readily calculate various quantities. The 
value p of the ground state phase is determined from 
7gl(V') = C^glCV') = 0, he., 






(All) 


The self-generated flux in the ground state is given by 
Eq. (19) and for o —>■ Oc can be approximated by 


27r 


£gl 

$0 


/3e 


Pt. / -«1 ^ 

l+/3i V 26o 


(A12) 
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where we took into account that for small ip, according 
to Eq. (2), (pi Ri ip/{l + /3i) and (p 2 ^ '0/(1 + a/ 32 )- 
The barrier between two wells is given by 

= (A13) 

The depinning current 7 c- can be found as an ex¬ 
tremum of 7 ( 0 ). The extremum is reached at 0 = 0dep 
satisfying the equation 7'(0dep) = 2a -|- 1260^gp = 0 . 
Thus, 



The value of 7 c- = 7(0dep) is 



7- = 0dep 




(A14) 


(A15) 


The value of 7 c+ cannot be calculated in the framework 
of our GL-approximation as 7 c+ corresponds to a large 
depinning phase. 

The eigenfrequency in each of the wells can be 
calculated as 


Wo = U"{ip) = = —4a Ri —4aie. (A16) 
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